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Hyper-Kamiokande will be a next generation underground water Cherenkov detec¬ 
tor with a total (fiducial) mass of 0.99 (0.56) million metric tons, approximately 20 
(25) times larger than that of Super-Kamiokande. One of the main goals of Hyper- 
Kamiokande is the study of CP asymmetry in the lepton sector using accelerator 
neutrino and anti-neutrino beams. 

In this paper, the physics potential of a long baseline neutrino experiment using the 
Hyper-Kamiokande detector and a neutrino beam from the J-PARC proton synchrotron 
is presented. The analysis uses the framework and systematic uncertainties derived from 
the ongoing T2K experiment. With a total exposure of 7.5 MW x 10 7 sec integrated 
proton beam power (corresponding to 1.56 x 10 22 protons on target with a 30 GeV 
proton beam) to a 2.5-degree off-axis neutrino beam, it is expected that the leptonic CP 
phase Scp can be determined to better than 19 degrees for all possible values of 5cp, and 
CP violation can be established with a statistical significance of more than 3cr (5<r) for 
76% (58%) of the Scp parameter space. Using both v e appearance and v jL disappearance 
data, the expected la uncertainty of sin 2 d 2 3 is 0.015(0.006) for sin 2 (9 23 = 0.5(0.45). 
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1. Introduction 

The discovery of neutrino oscillations by the Super-Kamiokande (Super-K) experiment in 
1998 [1] opened a new window to explore physics beyond the Standard Model (BSM). Evi¬ 
dence of neutrino oscillations is one of the most convincing experimental proofs known 
today for the existence of BSM physics at work. The mixing parameters of neutrinos were 
found to be remarkably different from those of quarks, which suggests the presence of an 
unknown flavor symmetry waiting to be explored. The extremely small masses of neutrinos 
compared with those of their charged partners lead to the preferred scenario of a seesaw 
mechanism [2-5], in which small neutrino masses are a reflection of the ultra-high energy 
scale of BSM physics. 

Furthermore, a theoretical framework called leptogenesis points to the intriguing possibility 
that CP asymmetries related to flavor mixing among the three generations of neutrinos may 
have played an important role in creating the observed matter-antimatter asymmetry in the 
universe [ 6 ]. This makes a study of the full picture of neutrino masses and mixings and the 
measurement of the CP asymmetry in the neutrino sector among the most important and 
urgent subjects in today’s elementary particle physics world. 

CP asymmetry in the neutrino sector arising from the presence of the phase which corre¬ 
sponds to the Kobayashi-Maskawa phase [7] in the quark sector, can only be seen if all the 
three mixing angles governing neutrino oscillations differ from zero. The Super-K detector 
has successfully measured all three angles. The angle 823 was first measured in atmospheric 
neutrino observations [ 1 ], 9 12 was constrained in solar neutrino observations [8] (together 
with another water Cherenkov detector SNO [9]), and the evidence of non-zero #13 was 
found by T2K [10] which used Super-K as the far detector. In 2013, T2K established /.y, —»• v e 
oscillation with 7.3 a significance, leading the way towards CP violation measurements in 
neutrinos [ 11 ] in combination with precise measurements of 9 13 by reactor neutrino exper¬ 
iments [12-14]. The highly successful Super-K program indicates that Hyper-Kamiokande 
(Hyper-K) is well placed to discover CP violation. 

In this paper, the physics potential of a long baseline neutrino experiment using the 
Hyper-Kamiokande detector and a neutrino beam from the J-PARC proton synchrotron 
is presented. 

The Hyper-K detector is designed as a next generation underground water Cherenkov 
detector that serves as a far detector of a long baseline neutrino oscillation experiment for the 
J-PARC neutrino beam and as a detector capable of observing proton decays, atmospheric 
and solar neutrinos, and neutrinos from other astrophysical origins. The baseline design of 
Hyper-K is based on the well-proven technologies employed and tested at Super-K. Hyper-K 
consists of two cylindrical tanks lying side-by-side, the outer dimensions of each tank being 
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Fig. 1: Schematic view of the Hyper-Kamiokande detector. 


48 (W) x 54 (H) x 250 (L) m 3 . The total (fiducial) mass of the detector is 0.99 (0.56) million 
metric tons, which is about 20 (25) times larger than that of Super-K. A proposed location 
for Hyper-K is about 8 km south of Super-K (and 295 km away from J-PARC) and 1,750 
meters water equivalent (or 648 m of rock) deep. The inner detector region is viewed by 
99,000 20-inch PMTs, corresponding to the PMT density of 20% photo-cathode coverage 
(the same as the second phase of Super-K). The schematic view of the Hyper-K detector is 
illustrated in Fig. 1. 

In addition to the long baseline neutrino oscillation experiment that is the main focus of 
this paper, Hyper-K will provide a rich program in a wide range of science [15]. The scope 
of the project includes observation of atmospheric and solar neutrinos, proton decays, and 
neutrinos from other astrophysical origins. The physics potential of Hyper-K is summarized 
in Table 1. 

2. Neutrino Oscillations and CP Violation 

2.1. Neutrino Oscillations in Three Flavor Framework 

Throughout this paper, unless stated otherwise, we consider the standard three flavor 
neutrino framework. The 3x3 unitary matrix U which describes the mixing of neutri¬ 
nos [16] (that is often referred to as the Maki-Nakagawa-Sakata-Pontecorvo (MNSP) or 
Maki-Nakagawa-Sakata (MNS) [16,17] matrix) relates the flavor and mass eigenstates of 
neutrinos: 

3 

l'a. = T: UgiVi- (a = e, Hi t), (1) 

i=l 

where n a (a = e,/r, r) and vfy = 1,2,3) denote, respectively, flavor and mass eigenstates of 
neutrinos. Using the standard parametrization, which can be found, e.g. in Ref. [18], U can 
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Table 1: Physics targets and expected sensitivities of the Hyper-Kamiokande experiment, based on the study shown in [15] except for the 
long baseline experiment that is described in this paper. Improvement is expected with further optimization of the detector design and 
development of reconstruction/analysis tools. Also, only selected values are listed; for example, other channels will be accessible for nucleon 
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be expressed as 
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where C{j = cos Oij, Sij = sin 0 t] , and 5cp — often called the Dirac CP phase —, is the 
Kobayashi-Maskawa type CP phase [7] in the lepton sector. On the other hand, the two 
phases, a 21 and 031 , — often called Majorana CP phases — exist only if neutrinos are of 
Majorana type [19-21], While the Majorana CP phases can not be observed in neutrino 
oscillation, they can be probed by lepton number violating processes such as neutrinoless 
double beta decay. 

In vacuum, the oscillation probability of u a —> up (a, f3 = e, fj,, r) for ultrarelativistic 
neutrinos is given by, 


P{u a -* up) 






i= 1 


2 


= 5 a p-4 J2 ^(K.UajUpiU^) sin 2 

i>j 



+2Y J ^{KiU aj Up i U*p j ) S \n 

i>j 



( 3 ) 


where E u is the neutrino energy, L is the baseline, Am| = rnf — m 2 (■ i,j = 1, 2,3) are the 
mass squared differences with m,i and rrij being the neutrino masses. For the CP conjugate 
channel, u Q — > vp, the same expression in Eq. (3) holds, but the matrix U is replaced by its 
complex conjugate (or equivalently 5cp — > —£>cp in Eq. (2)), resulting in the third term in 
Equation 3 switching sign. For neutrinos traveling inside matter, coherent forward scattering 
induces an asymmetry between the oscillation probabilities of neutrinos and antineutrinos 
supplementary to the intrinsic CP violation. 

The magnitude of the CP violation in neutrino oscillation can be characterized by the 
difference of probabilities between neutrino and anti-neutrino channels, which, in vacuum, 
is given by [22,23], 


AP q/3 = P{v a ->• up) - P(u a ->• up) = IQ J a p sin A 2 i sin A 32 sin A 3i , (4) 


and 


Jct/3 = A(7/q 1 [7*2 Up [ Up 2 ) = ±JcP, JcP — Sl2Cl2S23C23«13Ci 3 Sm ScP (5) 

with positive (negative) sign for (anti-)cyclic permutation of the flavor indices e, /j and r. 
The parameter Jqp is the lepton analogue of the CP-invariant factor for quarks, the unique 
and phase-convention-independent measure for CP violation [24]. Using the current best 
fitted values of mixing parameters [25], we get Jcp — 0.034 sin 6cp, or 

A P a p ~ ±0.55 sin <5 cp sin A 2 i sin A 3 2 sin A 3 i. ( 6 ) 

Thus, a large CP violation effects are possible in the neutrino oscillation. 
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In general, it is considered that CP violation in the neutrino sector which can be observed 
in the low energy regime, namely, in neutrino oscillation, does not directly imply the CP 
violation required at high energy for the successful leptogenesis in the early universe. It 
has been discussed, however, that they could be related to each other and the CP violating 
phase in the MNS matrix could be responsible also for the generation of the observed baryon 
asymmetry through leptogenesis in some scenarios. For example, in [26,27], in the context of 
the seesaw mechanism, it has been pointed out that assuming the hierarchical mass spectrum 
for right handed Majorana neutrinos with the lightest mass to be < 5 x 10 12 GeV, observed 
baryon asymmetry could be generated through the leptogenesis if | sin $13 sin^cpl >0.1, 
which is compatible with the current neutrino data. Hence, measurement of CP asymmetry 
in neutrino oscillations may provide a clue for understanding the origin of matter-antimatter 
asymmetry of the Universe. 

Since there are only three neutrinos, only two mass squared differences, and Am^, 

for example, are independent. Therefore, for a given energy and baseline, there are six inde¬ 
pendent parameters, namely, three mixing angles, one CP phase, and two mass squared 
differences, in order to describe neutrino oscillations. Among these six parameters, 612 and 
Ahave been measured by solar [9, 28, 29] and reactor [30-32] neutrino experiments. 
The parameters 622 , and | Am] 2 1 (only its absolute value) have been measured by atmo¬ 
spheric [33,34] and accelerator [35-38] neutrino experiments. Reactor experiment also starts 
to measure the atmospheric mass squared difference, [AmJjJ though the uncertainty is still 
larger [12]. Recently, $13 has also been measured by accelerator [10,11,39,40] and reactor 
experiments [12-14,41,42], The relatively large value of 6U3 opens the window to explore the 
CP phase (5cp) and the mass hierarchy (the sign of Am^) using neutrino oscillation. 

2.2. Physics Case with —» v e Oscillation 

The oscillation probability from to v e in accelerator experiments is expressed, to the first 
order of the matter effect, as follows [43]: 

P(V/i Ve) = 4c 2 3 S 2 3 S23 ' S™ 2 A 31 

+ 8 c 2 3 si 2 si 3 S 23 (ci 2 C 23 CosAcp - S12S13.S23) ■ cos A 32 • sinA 3 i • sin A21 
- 8 c 2 3 ci 2 C 23 Si 2 Si 3 'S 23 sinc>cp • sin A32 • sinA 3 i • sin A 2 i 
+ 4 6’ 2 2CJ3( C 2 2C23 + •Sl 2 s 23 s 13 — 2 c 12 C 23 s 12 S 23 s 13 COS 5 cp) ' sill 2 A21 
-84*13*23 • ^-(1 - 24s) • cos A32 • sinA 3 i 

+8Ci3Si34^ m 2 (1 _ 2 4) ■ s i n " A 31 > (7) 

where Ay is Am| L/4E/, and a = 2\Z2GFn e E u = 7.56 x 10 _5 [eV 2 ] x p[g/cm 3 ] x E u [ GeV]. 
The corresponding probability for a —> V e transition is obtained by replacing Scp —>• — 5cp 
and a —> —a. The third term, containing sin 5cp, is the CP violating term which flips sign 
between u and v and thus introduces CP asymmetry if sinAcp is non-zero. The last two 
terms are due to the matter effect. Those terms which contain a change their sign depending 
on the mass hierarchy. As seen from the definition of a, the amount of asymmetry due to the 
matter effect is proportional to the neutrino energy at a fixed value of L / E v . A direct test 
of CP violation, in a model independent way, is possible by measuring both neutrino and 
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antineutrino appearance probabilities. If the mass hierarchy is not known, the sensitivity of 
CP violation is affected by the presence of the matter effect. However, the mass hierarchy 
could be determined by the atmospheric neutrino measurement in Hyper-K and several 
measurements by other experiments. 

Currently measured value of 023 is consistent with maximal mixing, 023 ~ 7 t /4 [38,44,45]. 
It is of great interest to determine if sin 2 2023 is maximal or not, and if not 023 is less or 
greater than 7 r/ 4 , as it could constrain models of neutrino mass generation [46-51]. When 
we measure 023 with the survival probability P{y^ —)• v^ L ) which is proportional to sin 2 2023 
to first order, 


P(v,i -> *V) - 1 - 4c? 3 sl 3 [l - c 2 13 s 2 23 ] sin 2 (A m\ 2 L/AE V ) 
~ 1 — sin 2 2023 sin 2 (A?n 3 2 L/4E U ), 


(for ci 3 ~ 1 ) 


( 8 ) 

(9) 


there is an octant ambiguity: either 023 < 45° (in the first octant) or 6-23 > 45° (in the second 
octant). By combining the measurement of P(z/^ —> v e ), the 623 octant can be determined. 

2.3. Anticipated Neutrino Physics Landscape in the 2020s and Uniqueness of This 
Experiment 

Before Hyper-K commences data taking in ~ 2025, we expect a number of ongoing and 
planned neutrino experiments as well as cosmological observations will advance our under¬ 
standing of neutrino physics. In addition to accelerator and reactor experiments, Super-K 
will provide precise measurements of neutrino oscillation parameters from atmospheric neu¬ 
trino observations, and will look for the mass hierarchy and the octant of 023. Cosmological 
observations will provide the information on neutrino masses. An observation of neutrino¬ 
less double /? decay in the next 10 years would be evidence that the neutrino is a Majorana 
particle with the inverted mass hierarchy. Following this progress, we definitely need a new 
experiment to discover CP violation in neutrinos, and to unambiguously establish the mass 
hierarchy and 623 octant. For these purposes, we propose the Hyper-K experiment with the 
J-PARC neutrino beam. 

The Hyper-K experiment will have several unique advantages. 

o The experiment will have high statistics of neutrino events thanks to the large fiducial 
mass and the high power J-PARC neutrino beam, 
o The relatively short baseline among the proposed long baseline experiments results in a 
small ambiguity from the matter effect. 

o The experiment will operate in the same beam line as T2K with the same off-axis 
configuration. The features of the neutrino beam and the operation of the high power 
beam are well understood. 

o The systematc errors are already well understood based on Super-K and T2K, allowing 
reliable extrapolations. 

With these features, Hyper-K will be one of the most sensitive experiments to probe neutrino 
CP violation, as we present in this paper. 
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Table 2: Planned parameters of the J-PARC Main Ring for fast extraction. Numbers in 
parentheses are those achieved up until May 2013. 


Parameter 


Value 

Circumference (m) 


1567.5 

Kinetic energy (GeV) 


30 

Beam intensity (ppp) 

2.0 x 10 

14 (1.24 x 10 14 ) 

(ppb) 

2.5 x 10 

13 (1.57 x 10 13 ) 

Harmonic number 


9 

Number of bunches per spill 


8 

Spill width (/j, s) 


~ 5 

Bunch full width at extraction (ns) 


~ 50 

Maximum RF voltage (kV) 

560 

(280) 

Repetition period (sec) 

1.28 

(2.48) 

Beam power (kW) 

750 

(240) 


3. Experimental setup 

3.1. J-PARC accelerator and neutrino beamline 

An intense and high quality neutrino beam is a key for the success of a long baseline neu¬ 
trino oscillation experiment. J-PARC (Japan Proton Accelerator Research Complex) is one 
of world leading facilities in neutrino physics, currently providing a beam for the T2K exper¬ 
iment. We will utilize the full potential of this existing facility with future increase of the 
beam power to the design value of 750 kW and beyond. 

The J-PARC accelerator cascade [52] consists of a normal-conducting LINAC as an injec¬ 
tion system, a Rapid Cycling Synchrotron (RCS), and a Main Ring synchrotron (MR). In 
the fast extraction mode operation, MR has achieved 1.24xl0 14 protons per pulse (ppp) 
beam intensity, which is a world record for extracted ppp for any synchrotron. The corre¬ 
sponding beam power is 240 kW. The upgrade scenario of J-PARC accelerator [53] is being 
implemented to reach the design power of 750 kW in forthcoming years, with a typical 
planned parameter set as listed in Table 2. This will double the current repetition rate by (i) 
replacing the magnet power supplies, (ii) replacing the RF system, and (iii) upgrading injec¬ 
tion/extraction devices. The design power of 750 kW will be achieved well before Hyper-K 
will start data taking. Furthermore, conceptual studies on how to realize 1~2MW beam 
powers and even beyond are now underway [54], such as by raising the RCS top energy, 
enlarging the MR aperture, or inserting an “emittance-damping” ring between the RCS and 
MR. 

Figure 2 shows an overview of the neutrino experimental facility [55,56]. The primary 
beamline guides the extracted proton beam to a production target/pion-focusing horn system 
in a target station. The pions decay into muons and neutrinos during their flight in a 110 
m-long decay volume. A graphite beam dump is installed at the end of the decay volume, 
and muon monitors downstream of the beam dump monitor the muon profile. A neutrino 
near detector complex is situated 280 m downstream of the target to monitor neutrinos 
at production. To generate a narrow band neutrino beam, the beamline utilizes an off-axis 
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Fig. 2: The neutrino experimental facility (neutrino beamline) at J-PARC. 


Table 3: Acceptable beam power and achievable parameters for each beamline compo¬ 
nent [58]. Limitations as of May 2013 are also given in parentheses. 


Component 

Beam power/parameter 

Target 

3.3X10 14 ppp 

Beam window 

3.3xl0 14 ppp 

Horn 



Cooling for conductors 


2 MW 

Stripline cooling 

1~2 MW 

( 400 kW ) 

Hydrogen production 

1~2 MW 

( 300 kW ) 

Horn current 

320 kA 

( 250 kA ) 

Power supply repetition 

1 Hz 

( 0.4 Hz ) 

Decay volume 


4 MW 

Hadron absorber/beam dump 


3 MW 

Water cooling facilities 

~2 MW 

( 750 kW ) 

Radiation shielding 

4 MW 

( 750 kW ) 

Radioactive air leakage to the target station ground floor 

~2 MW 

( 500 kW ) 

Radioactive cooling water treatment 

~2 MW 

( 600 kW ) 


beam configuration [57], with the capability to vary the off-axis angle in the range from 2.0° 
to 2.5°. The latter value has been used for the T2K experiment and is assumed also for the 
proposed project. The centerline of the beamline extends 295 km to the west, passing midway 
between Tochibora (Hyper-K candidate site) and Mozumi (where Super-K is located), so 
that both sites have identical off-axis angles. 

Based on the considerable experience gained on the path to achieving 240 kW beam power 
operation, improvement plans to realize 750 kW operation, such as improving the activated 
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air confinement in the target station and expanding the facilities for the treatment of acti¬ 
vated water, are being implemented and/or proposed. Table 3 gives a summary of acceptable 
beam power and/or achievable parameters for each beamline component [58], after the 
proposed improvements in forthcoming years. 


3.2. Near detectors 

The accelerator neutrino event rate observed at Hyper-K depends on the oscillation proba¬ 
bility, neutrino flux, neutrino interaction cross-section, detection efficiency, and the detector 
fiducial mass of Hyper-K. To extract estimates of the oscillation parameters from data, one 
must model the neutrino flux, cross-section and detection efficiency with sufficient precision. 
In the case of the neutrino cross-section, the model must describe the exclusive differential 
cross-section that includes the dependence on the incident neutrino energy, E v , the kinemat¬ 
ics of the outgoing lepton, pi and 6 i, and the kinematics of final state hadrons and photons. 
In our case, the neutrino energy is inferred from the lepton kinematics, while the modeling 
of reconstruction efficiencies depends on the hadronic final state as well. 

The neutrino flux and cross-section models can be constrained by data collected at near 
detectors, situated close enough to the neutrino production point so that oscillation effects 
are negligible. Our approach to using near detector data will build on the experience of 
T2K while considering new near detectors that may address important uncertainties in the 
neutrino flux or cross-section modeling. 

The conceptual design of the near detectors is being developed based on the physics sensi¬ 
tivity studies described in Section 4. In this section, we present basic considerations on the 
near detector requirements and conceptual designs. More concrete requirements and detector 
design will be presented in future. 

We assume to use T2K near detectors [55], INGRID and ND280, possibly with an upgrade. 
The INGRID detector [59] consists of 16 iron-scintillator modules configured in a cross 
pattern centered on the beam axis 280 m downstream from the T2K target. The rate of 
interactions in each module is measured and a profile is constructed to constrain the neu¬ 
trino beam direction. The ND280 off-axis detector is located 280 m downstream from the 
T2K target as well, but at an angle of 2.5 degrees away from the beam direction. The POD 
7r° detector [60], time projection chambers (TPCs) [61], fine grain scintillator bar detectors 
(FGDs) [62] and surrounding electromagnetic calorimeters (ECALs) [63]. The detectors are 
immersed in a 0.2 T magnetic held and the magnetic yoke is instrumented with plastic 
scintillator panels for muon range detection [64], The magnetic held allows for momentum 
measurement and sign selection of charged particles. The magnetization of ND280 is par¬ 
ticularly important for operation in antineutrino mode where the neutrino background is 
large. In that case, ND280 is able to separate the “right-sign” p + from the “wrong-sign” 

. The POD and FGDs act as the neutrino targets, while the TPCs provide measurements 
of momentum and ionizing energy loss for particle identification. The POD and one of the 
FGDs include passive water layers that allow for neutrino interaction rate measurements 
on the same target as Super-K. ND280 has been employed to measure the rates of charged 
current z/ /( and u e interactions, as well as NC7 t° interactions. 

The T2K collaboration is in the process of discussing various upgrade possibilities at the 
ND280 site [65]. These include the deployment of heavy water (D2O) within the passive 
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water targets in FGD2 that would allow the extraction of neutrino interaction properties 
on the quasi-free neutron in deuterium via a subtraction with data taken with light water 
H2O. The use of a water-based liquid scintillator (WbLS) developed at BNL [66] is being 
explored in the context of a tracking detector with comparable or finer granularity than the 
FGD to allow the detailed reconstruction of hadronic system emerging from the neutrino 
interactions or a larger detector with coarser segmentation that would allow high statistics 
studies. Either would significantly enhance the study of neutrino interactions on water by 
reducing the reliance on subtraction and enhancing the reconstruction capabilities relative 
to the currently deployed passive targets. Finally, a high pressure TPC that can contain 
various noble gases (He, Ne, Ar) to serve both as the target and tracking medium is being 
studied. Such a detector would allow the ultimate resolution of the particles emitted from 
the target nucleus while allowing a study of the A-dependence of the cross-sections and final 
state interactions to rigorously test models employed in neutrino event generators. 

Since many of the uncertainties on the modeling of neutrino interactions arise from uncer¬ 
tainties on nuclear effects, the ideal near detector should include the same nuclear targets as 
the far detector. In T2K near detectors, the POD [60] and FGD [62] detectors include passive 
water layers, however extracting water only cross sections requires complicated analyses that 
subtract out the interactions on other materials in the detectors. An alternative approach 
is to build a water Cherenkov (WC) near detector to measure the cross section on H2O 
directly and with no need for a subtraction analysis. This approach was taken by K2K [35] 
and was proposed for T2K [67]. The MiniBooNE experiment has also employed a mineral 
oil Cherenkov detector at a short baseline to great success [68]. A WC near detector design 
is largely guided by two requirements: 

(1) The detector should be large enough to contain muons up to the momentum of interest 
for measurements at the far detector, and to provide sufficient radiation length for 
detection of gamma rays. 

(2) The detector should be far enough from the neutrino production point so that there 
is minimal pile-up of interactions in the same beam timing bunch. 

These requirements lead to designs for kiloton size detectors located at intermediate 
distances, 1-2 km from the target, for the J-PARC neutrino beam. 

The main disadvantage of the WC detector is the inability to separate positively and 
negatively charged leptons, and hence antineutrino and neutrino interactions. This ability 
is especially important for a CP asymmetry measurement where the wrong sign contri¬ 
bution to the neutrino flux should be well understood. Hence, the WC detector will most 
likely be used in conjunction with a magnetized tracking detector such as ND280. Recent 
developments in the addition of Gadolinium (Gd) [69] and Water-based Liquid Scintilla¬ 
tor (WbLS) compounds [66] to water do raise the possibility to separate neutrino and 
antineutrino interactions by detecting the presence of neutrons or protons in the final state. 

Two conceptual designs for possible intermediate WC detectors have been studied. Unoscil¬ 
lated Spectrum (TITUS) is a 2 kiloton WC detector located about 2 km from the target 
at the same off-axis angle as the far detector. At this baseline the detector sees fluxes for 
the neutral current and u e backgrounds that are nearly identical to the Hyper-K fluxes. The 
detector geometry and the presence of a muon range detector are optimized to detect the 
high momentum tail of the muon spectrum. The use of Gd in TITUS to separate neutrino 
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Table 4: Parameters of the Hyper-Kamiokande baseline design. 


Detector type 


Ring-imaging 

water Cherenkov detector 

Candidate site 

Address 

Lat. 

Long. 

Alt. 

Overburden 

Cosmic Ray Muon flux 

Off-axis angle for the J-PARC u 
Distance from the J-PARC 

Tochibora mine 

Kamioka town, Gifu, JAPAN 
36°21 / 20.105"N t 
137°18'49.137"E t 

508 m 

648 m rock 

(1,750 m water equivalent) 

~ 8 X 10 - ' sec -1 cnK 2 

2.5° (same as Super-K) 

295 km (same as Super-K) 

Detector geometry 

Total Water Mass 

Inner Detector (Fiducial) Mass 
Outer Detector Mass 

0.99 Megaton 

0.74 (0.56) Megaton 

0.2 Megaton 

Photo-sensors 

Inner detector 

Outer detector 

99,000 20-inch cp PMTs 

20 % photo-coverage 

25,000 8-inch cp PMTs 

Water quality 

light attenuation length 

Rn concentration 

> 100 m @ 400 nm 
< 1 lnBq/m 3 


' World geographical coordination system 


and antineutrino interactions is being studied. The zzPRISM detector is located 1 km from 
the target and is 50 m tall, covering a range of off-axis angles from 1-4 degrees. The z'PRISM 
detector sees a range of neutrino spectra, peaked at energies from 0.4 to 1.0 GeV depending 
on the off-axis angle. The purpose of zaPRISM is to use these spectra to better probe the 
relationship between the incident neutrino energy and final state lepton kinematics, a part 
of the interaction model with larger uncertainties arising from nuclear effects. 

3.3. Hyper-Kamiokande 

Hyper-Kamiokande is to be the third generation water Cherenkov detector in Kamioka, 
designed for a wide variety of neutrino studies and nucleon decay searches. Its total (fiducial) 
water mass of one (0.56) million tons would be approximately 20 (25) times larger than that 
of Super-Kamiokande. Table 4 summarizes the baseline design parameters of the Hyper-K 
detector. 

In the baseline design, the Hyper-K detector is composed of two separated caverns as 
shown in Fig. 1, each having an egg-shape cross section 48 meters wide, 54 meters tall, and 
250 meters long as shown in Fig. 3 and 4. The welded polyethylene tanks are filled up to a 
depth of 48 m with ultra-pure water: the total water mass equals 0.99 million tons. 

Each tank will be optically separated by segmentation walls located every 49.5 m to form 
5 (in total 10) compartments as shown in Fig. 4, such that event triggering and event 
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CROSS SECTION 


Measurement Facility Area 



Fig. 3: Cross section view of the Hyper-Kamiokande detector. 


CONCEPT 



5«49500(Pitch of Interamedatn Cross Walls)=24750O 


Fig. 4: Profile of the Hyper-K detector. Top: the detector segmentation. Bottom: PMT arrays 
and the support structure for the inner and outer detectors. Each quasi-cylindrical tank lying 
horizontally is segmented by intermediate walls into five compartments. 
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reconstruction can be performed in each compartment separately and independently. Because 
the compartment dimension of 50 m is comparable with that of Super-K (36 m) and is shorter 
than the typical light attenuation length in water achieved by the Super-K water filtration 
system (> 100 m @ 400 nm), we expect that the detector performance of Hyper-K for beam 
and atmospheric neutrinos will be effectively the same as that of Super-K. 

The water in each compartment is further optically separated into three regions. The inner 
region has a barrel shape of 42 m in height and width, and 48.5 m in length, and is viewed by 
an inward-facing array of 20-inch diameter photomultiplier tubes (PMTs). The entire array 
consists of 99,000 Hamamatsu R3600 PMTs, uniformly surrounding the region and giving 
a photocathode coverage of 20%. The PMT type, size, and number density are subject to 
optimization. We have been also developing new photosensors as possible alternative options 
to the R3600, such as a PMT with a box-and-line dynode and a hybrid photo-detector (HPD), 
both with a high quantum efficiency photocathode. An outer region completely surrounds 
the 5 (in total 10) inner regions and is equipped with 25,000 8-inch diameter PMTs. This 
region is 2 m thick at the top, bottom, and barrel sides, except at both ends of each cavern, 
where the outer region is larger than 2 nr due to rock engineering considerations. A primary 
function of the outer detector is to reject entering cosnric-ray muon backgrounds and to help 
in identifying nucleon decays and neutrino interactions occurring in the inner detector. The 
middle region or dead space is an uninstrumented, 0.9 nr thick shell between the inner and 
outer detector volumes where the stainless steel PMT support structure is located. Borders 
of both inner and outer regions are lined with opaque sheets. This dead space, along with 
the outer region, acts as a shield against radioactivity from the surrounding rock. The total 
water mass of the inner region is 0.74 million tons and the total fiducial mass is 10 times 
0.056 = 0.56 million tons. The fiducial volume is defined as the region formed by a virtual 
boundary located 2 m away from the inner PMT plane. 

The estimated cosmic-ray muon rate around the Hyper-K detector candidate site is ~ 8 
x 10 -7 sec _1 cm~ 2 which is roughly 5 times larger than the flux at Super-K’s location (~ 
1.5 x ICE 7 sec _1 cm -2 ). The expected deadtime due to these muons is less than 1% and 
negligible for long baseline experiments, as well as nucleon decay searches and atmospheric 
neutrino studies. 

Water is the target material and signal-sensitive medium of the detector, and thus its 
quality directly affects the physics sensitivity. In Super-Kamiokande the water purification 
system has been continually modified and improved over the course of two decades. As a 
result, the transparency is now kept above 100 m and is very stable, and the radon concen¬ 
tration in the tank is held below 1 mBq/m 3 . Following this success, the Hyper-Kamiokande 
water system has been designed based on the current Super-Kamiokande water system with 
scaling up the process speeds to 1200 m 3 /hour for water circulation and 400 m 3 /hour for 
radon free air generation. With these systems, the water quality in Hyper-Kamiokande is 
expected to be same as that in Super-Kamiokande. Adding dissolved gadolinium sulfate 
for efficient tagging of neutrons has been studied as an option to enhance Hyper-K physics 
capability. The feasibility of adding Gd to Super-K [69] is now under study with EGADS 
(Evaluating Gadolinium’s Action on Detector Systems) project in Kamioka. We have been 
careful to keep the possibility of gadolinium loading in mind when designing the overall 
Hyper-Kamiokande water system. 
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Fig. 5: PID likelihood functions for electron (blue solid histogram) and // (red dashed his¬ 
togram) with 500MeV/c momentum. A negative (positive) value indicates electron-like 
(//-like) particle. 


We have evaluated the expected performance of the Hyper-K detector using the MC sim¬ 
ulation and reconstruction tools under development. We have been developing a detector 
simulation dedicated to Hyper-K based on “WCSim,” [70] which is an open-source water 
Cherenkov detector simulator based on the GEANT4 library [71,72]. A new reconstruc¬ 
tion algorithm developed for Super-K/T2K [11], named “fiTQun,” has been adopted for the 
Hyper-K analysis. It uses a maximum likelihood fit with charge and time probability density 
functions constructed for every PMT hit assuming several sets of physics variables (such as 
vertex, direction, momentum, and particle type) [11,73]. 

As an example of the evaluation, electrons and muons with 500 MeV/c are generated with 
a fixed vertex (at the center of the tank) and direction (toward the barrel of the tank) in the 
Hyper-K detector simulation. Figure 5 shows the likelihood function for the particle identi¬ 
fication. A negative (positive) value indicates electron-like (//-like) particle. It demonstrates 
a clear separation of electrons and muons. The obtained performance of Hyper-Kamiokande 
is compared with the performance of SK-II (20% photo coverage, old electronics) and SK-IV 
(40% photo coverage, new electronics) in Table 5. The vertex resolution for muon events will 
be improved to the same level as Super-K with an update of the reconstruction program. 
From the preliminary studies, the performance of Hyper-K is similar to or possibly better 
than SK-II or SK-IV with the new algorithm. In the physics sensitivity study described 
in Section 4, a Super-K full MC simulation with the SK-IV configuration is used because 
it includes the simulation of new electronics and is tuned with the real data, while giving 
similar performance with Hyper-K as demonstrated above. 

4. Physics Sensitivities 

4-1. Overview 

As discussed in Sec. 2 . 2 , a comparison of muon-type to electron-type transition probabilities 
between neutrinos and anti-neutrinos is one of the most promising methods to observe the 
lepton CP asymmetry. Recent observation of a nonzero, rather large value of $13 [10,14,41,42] 
makes this exciting possibility more realistic. 
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Table 5: Comparison of performance of SK-II (20% photo-coverage), SK-IV (40% photo¬ 
coverage), and the expected performance of Hyper-Kamiokande baseline design (20% photo¬ 
coverage) with preliminary Hyper-K simulation and reconstruction. 



SK-II 

SK-IV 

Hyper-K 

Particle type (p =500MeV/c) 

e 

l 1 

e 


e 

l 1 

Vertex resolution 

28 cm 

23 cm 

25 cm 

17 cm 

27 cm 

30 cm 

Particle identification 

98.5% 

99.0% 

98.8% 

99.5% 

>99.9% 

99.2% 

Momentum resolution 

5.6% 

3.6% 

4.4% 

2.3% 

4.0% 

2.6% 


neutrino 



anti-neutrino 



Fig. 6: Oscillation probabilities as a function of the neutrino energy for —> v e (left) and 

—> V e (right) transitions with L=295 km and sin 2 2013 = 0.1. Black, red, green, and blue 
lines correspond to dcp = 0, and — ^vr, respectively. Solid (dashed) line represents the 

case for a normal (inverted) mass hierarchy. 


Figure 6 shows the —>• v e and —> V e oscillation probabilities as a function of the true 
neutrino energy for a baseline of 295 km. The Earth matter density of 2.65/cm 3 is used in 
this analysis. The cases for Scp = 0, ^7r, it, and — ^vr, are overlaid. Also shown are the case of 
normal mass hierarchy (Am| 2 > 0) with solid lines and inverted mass hierarchy (Am| 2 < 0) 
with dashed lines. The oscillation probabilities depend on the value of <5 cp, and by comparing 
the neutrinos and anti-neutrinos, one can see the effect of CP violation. 

There are sets of different mass hierarchy and values of 5 cp which give similar oscillation 
probabilities. This is known as the degeneracy due to unknown mass hierarchy and may 
introduce an ambiguity if we do not know the true mass hierarchy. Because there are a 
number of experiments planned to determine mass hierarchy in the near future, it is expected 
that the mass hierarchy will be determined by the time Hyper-K starts to take data. If 
not, Hyper-K itself has a sensitivity to the mass hierarchy by the atmospheric neutrino 
measurements as shown in Table 1. Furthermore, a combined analysis of the accelerator and 
atmospheric neutrino data in Hyper-K will enhance the sensitivity as shown in Sec. 4.7. 
Thus, the mass hierarchy is assumed to be known in this analysis, unless otherwise stated. 

Figure 7 shows the contribution from each term of the u^ —>• v e oscillation probability for¬ 
mula, Eq.(7), for L = 295 km, sin 2 26*13 = 0.1, sin 2 26*23 = TO, Scp = tt/ 2, and normal mass 


18/40 
































0 1 2 
E v (GeV) 


Fig. 7: Oscillation probability of —> v e as a function of the neutrino energy with a baseline 
of 295 km. sin 2 2#i3 = 0.1, Scp = and normal hierarchy are assumed. Contribution from 
each term of the oscillation probability formula is shown separately. 


hierarchy. For E u ~ 0.6 GeV which gives sinA32 — sin A31 ~ 1, 

Pjvy. ->• Ve) - P{p» Ve) _ -16J CP sill A 2i + - 2 S j 3 ) 

P{Vp t r'e) 4“ P{y^i t ^e) ®*T3®13®23 

= —0.28 sin S + 0.07. (11) 

The effect of CP violating term can be as large as 28%, while the matter effect is much 
smaller. 

The uncertainty of Earth density between Tokai and Kamioka is estimated to be at most 
6% [74], Because the matter effect contribution to the total appearance probability is less 
than 10% for 295km baseline, the uncertainty from matter density is estimated to be less 
than 0.6% and neglected in this analysis. 

Due to the relatively short baseline and thus lower neutrino energy at the oscillation max¬ 
imum, the contribution of the matter effect is smaller for the J-PARC to Hyper-Kamiokande 
experiment compared to other proposed experiments like LBNE in the United States [75] 
or LBNO in Europe [76]. Thus the CP asymmetry measurement with J-PARC to Hyper-K 
long baseline experiment has less uncertainty related to the matter effect, while other exper¬ 
iments with > 1000 km baseline have much better sensitivity to the mass hierarchy with 
accelerator neutrino beams 1 . The sensitivities for CP violation and mass hierarchy can be 
further enhanced by combining measurements with different baseline. 

The analysis method is based on a framework developed for the sensitivity study by T2K 
reported in [77]. A binned likelihood analysis based on the reconstructed neutrino energy 
distribution is performed using both v e (v e ) appearance and (v^) disappearance samples 
simultaneously. In addition to sin 2 2#i3 and Scp, sin 2 823 and Am| 2 are also included as free 
parameters in the fit. Table 6 shows the nominal oscillation parameters used in the study 
presented in this paper, and the treatment during the fitting. Systematic uncertainties are 
estimated based on the experience and prospects of the T2K experiment, and implemented 
as a covariance matrix which takes into account the correlation of uncertainties. 


1 Note that Hyper-K has sensitivity to the mass hierarchy using atmospheric neutrinos as shown 
in Table 1. 
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Table 6: Oscillation parameters used for the sensitivity analysis and treatment in the fitting. 
The nominal values are used for figures and numbers in this section, unless otherwise stated. 


Parameter 

Nominal value 

Treatment 

sin 2 2013 

0.10 

Fitted 

5cp 

0 

Fitted 

sin 2 023 

0.50 

Fitted 

<M 

(MOO 

<1 

2.4 x 10 -3 eV 2 

Fitted 

Mass hierarchy 

Normal or Inverted 

Fixed 

sin 2 2012 

0.8704 

Fixed 

Am|i 

7.6 x 10 -5 eV 2 

Fixed 


An integrated beam power of 7.5 MWxlO' sec is assumed in this study. It corresponds to 
1.56 x 10 22 protons on target with 30GeV J-PARC beam. We have studied the sensitivity 
to CP violation with various assumptions of neutrino mode and anti-neutrino mode beam 
running time ratio for both normal and inverted mass hierarchy cases. The dependence of 
the sensitivity on the ratio is found to be not significant between v\i 7=1:1 to 1:5. In this 
paper, v.V ratio is set to be 1:3 so that the expected number of events are approximately 
the same for neutrino and anti-neutrino modes. 

4-2. Neutrino flux 

The neutrino flux is estimated by T2K collaboration [78] by simulating the J-PARC 
neutrino beam line while tuning the modeling of hadronic interactions using data from 
NA61/SHINE [79,80] and other experiments measuring hadronic interactions on nuclei. To 
date, NA61/SHINE has provided measurements of pion and kaon production multiplicities 
for proton interactions on a 0.04 interaction length graphite target, as well as the inelastic 
cross section for protons on carbon. Since “thin” target data are used, the secondary inter¬ 
actions of hadrons inside and outside of the target are modeled using other data or scaling 
the NA61/SHINE data to different center of mass energies or target nuclei. NA61/SHINE 
also took data with a replica of 90cm-long T2K target, which will reduce the uncertainties 
related to the secondary interactions inside of the target. 

For the studies presented in this document, the T2K flux simulation has been used with 
the horn currents raised from 250 kA to 320 kA. The flux is estimated for both polarities of 
the horn fields, corresponding to neutrino enhanced and antineutrino enhanced fluxes. The 
calculated fluxes at Hyper-K, without oscillations, are shown in Fig. 8. 

The sources of uncertainty in the T2K flux calculation include: 

o Uncertainties on the primary production of pions and kaons in proton on carbon 
collisions. 

o Uncertainties on the secondary hadronic interactions of particles in the target or beam 
line materials after the initial hadronic scatter. 

o Uncertainties on the properties of the proton beam incident on the target, including the 
absolute current and the beam profile. 

o Uncertainties on the alignment of beam line components, including the target and 
magnetic horns. 
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Hyper-K Flux for Neutrino Mode Hyper-K Flux for Antieutrino Mode 



Fig. 8: The predicted Hyper-K neutrino fluxes from the J-PARC beam without oscillations. 
The neutrino enhanced beam is shown on the left and the antineutrino enhanced beam is 
shown on the right. 
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Fig. 9: The predicted uncertainty on the neutrino flux calculation assuming replica target 
hadron production data are available. 


o Uncertainties on the modeling of the horn fields, including the absolute field strength 
and asymmetries in the field. 

The uncertainties on the hadronic interaction modeling are the largest contribution to 
the flux uncertainty and may be reduced by using the hadron production data with a 
replica of T2K target. A preliminary analysis using a subset of the replica target data 
from NA61/SHINE has shown that it can be used to predict the T2K flux [81]. Since it is 
expected that replica target data will be available for future long baseline neutrino experi¬ 
ments, the Hyper-K flux uncertainty is estimated assuming the expected uncertainties on the 
measurement of particle multiplicities from the replica target. Hence, uncertainties related 
to the modeling of hadronic interactions inside the target are no longer relevant, however, 
uncertainties for interactions outside of the target are considered. The uncertainties on the 
measured replica target multiplicities are estimated by applying the same uncertainties that 
NA61/SHINE has reported for the thin target multiplicity measurements. 
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Fig. 10: The uncertainty on the far-to-near flux ratio for near detectors at 280 m, 1 km and 
2 km. Left: neutrino enhanced beam. Right: antineutrino enhanced beam. Top: the focused 
component of the beam. Bottom: the defocused component of the beam. 


The total uncertainties on the flux as function of the neutrino energy are shown in Fig. 9. 
In oscillation measurements, the predicted flux is used in combination with measurements 
of the neutrino interaction rate from near detectors. Hence, it is useful to consider the 
uncertainty on the ratio of the flux at the far and near detectors: 

' 5F ' m( - E '' , = 5 (St!) <i2) 

Here <f >and ^nd(^) are the predicted fluxes at Hyper-K and the near detector 
respectively. T2K uses the ND280 off-axis detector located 280 m from the T2K target. 
At that distance, the beam-line appears as a line source of neutrinos, compared to a point 
source seen by Hyper-K, and the far-to-near ratio is not flat. For near detectors placed 
further away, at 1 or 2 km for example, the far-to-near flux ratio becomes more flat and 
there is better cancellation of the flux uncertainties between the near and far detectors. 
Fig. 10 shows how the uncertainty on the far-to-near ratio evolves for baselines of 280 m, 
1 km and 2 km. While this extrapolation uncertainty is reduced for near detectors further 
from the production point, even the 280 m to Hyper-K uncertainty is less than 1% near the 
flux peak energy of 600 MeV. 
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Fig. 11: Reconstructed neutrino energy distribution of the v e candidate events. Normal mass 
hierarchy with sin 2 20 \s = 0.1 and Scp = 0 is assumed. 

4-3. Expected observables at Hyper-K 

Interactions of neutrinos in the Hyper-K detector are simulated with the NEUT program 
library [82-84], which is used in both Super-K and T2K. The response of the detector is 
simulated using the Super-K full Monte Carlo simulation based on the GEANT3 package [85]. 
The simulation is based on the SK-IV configuration with the upgraded electronics and DAQ 
system. Events are reconstructed with the Super-K reconstruction software. As described in 
Sec. 3.3, the performance of Hyper-K detector for neutrinos with J-PARC beam energy is 
expected to be similar to that of Super-K. Thus, the Super-K full simulation gives a realistic 
estimate of the Hyper-K performance. 

The criteria to select v e and candidate events are based on those developed for and 
established with the Super-K and T2K experiments. Fully contained (FC) events with a 
reconstructed vertex inside the fiducial volume (FV) and visible energy (E v ; s ) greater than 
30 MeV are selected as FCFV neutrino event candidates. In order to enhance charged current 
quasielastic (CCQE, vi + n —> l~ + p or Vi + p —> l + + n) interaction, a single Cherenkov ring 
is required. 

Assuming a CCQE interaction, the neutrino energy {E™ c ) is reconstructed from the energy 
of the final state charged lepton (Ed) and the angle between the neutrino beam and the 
charged lepton directions (9d) as 

rec = 2(m n - V)E e + m 2 - (m n - V) 2 - mj 
v 2 (m n -V - Eg + pe cos Qg) 

where m n ,m p ,mg are the mass of neutron, proton, and charged lepton, respectively, pg is 
the charged lepton momentum, and V is the nuclear potential energy (27MeV). 

Then, to select t'e/Te candidate events the following criteria are applied; 

o The reconstructed ring is identified as electron-like (e-like), 
o The visible energy (E v i s ) is greater than 100 MeV. 
o There is no decay electron associated to the event, 
o The reconstructed energy ( E ™ c ) is less than 1.25 GeV. 

o In order to reduce the background from mis-reconstructed 7r° events, additional criteria 
using a reconstruction algorithm recently developed for T2K (fiTQun, see Sec. 3.3) is 
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Table 7: The expected number of u e candidate events. Normal mass hierarchy with sin 2 26\^ = 
0.1 and 5cp = 0 are assumed. Background is categorized by the flavor before oscillation. 



signal 

Vft “S' Ve Vft -S' T e 

u„. CC 

1% CC 

BG 

v e CC F e CC 

NC 

BG Total 

Total 

v mode 

3016 

28 

11 

0 

503 

20 

172 

706 

3750 

v mode 

396 

2110 

4 

5 

222 

396 

265 

891 

3397 
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Fig. 12: Reconstructed neutrino energy distribution of the candidate events. 


Table 8: The expected number of zz /t candidate events. 



v tL CC 

CC 

u e CC 

v e CC 

NC 

v,i -S' V e 

total 

v mode 

17225 

1088 

11 

1 

999 

49 

19372 

v mode 

10066 

15597 

7 

7 

1281 

6 

26964 


applied. With a selection based on the reconstructed 7r° mass and the ratio of the best-fit 
likelihoods of the 7r° and electron fits as used in T2K [11], the remaining 7r° background 
is reduced to about 30% compared to the previous study [15]. 

Figure 11 shows the reconstructed neutrino energy distributions of v e events after all the 
selections. The expected number of v e candidate events is shown in Table 7 for each signal 
and background component. In the neutrino mode, the dominant background component 
is intrinsic v e contamination in the beam. The mis-identified neutral current 7T° production 
events are suppressed thanks to the improved 7r° rejection. In the anti-neutrino mode, in 
addition to v e and z7 /i; v e and components have non-negligible contributions due to larger 
fluxes and cross-sections compared to their counterparts in the neutrino mode. 

For the Vp/vp candidate events the following criteria are applied; 

o The reconstructed ring is identified as muon-like (/i-like). 

o The reconstructed muon momentum is greater than 200MeV/c. 

o There is at most one decay electron associated to the event. 
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Fig. 13: Top: Reconstructed neutrino energy distribution for several values of 5c p- sin 2 2(9 1,3 = 
0.1 and normal hierarchy is assumed. Bottom: Difference of the reconstructed neutrino 
energy distribution from the case with 5 cp = 0°. The error bars represent the statistical 
uncertainties of each bin. 


Neutrino mode: Disappearance Antineutrino mode: Disappearance 



Fig. 14: Reconstructed neutrino energy distribution of candidates for several values of 
5c P- 


Figure 12 shows the reconstructed neutrino energy distributions of the selected events. 

Table 8 shows the number of candidate events for each signal and background component. 
For the neutrino mode, most of the events are due to while in the anti-neutrino mode 
the contribution from wrong-sign components is significant. 

The reconstructed neutrino energy distributions of v e events for several values of 6 cp are 
shown in the top plots of Fig. 13. The effect of 5cp is clearly seen using the reconstructed 
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neutrino energy. The bottom plots show the difference of reconstructed energy spectrum from 
5cp = 0° f° r the cases 5 = 90°, —90° and 180°. The error bars correspond to the statistical 
uncertainty. By using not only the total number of events but also the reconstructed energy 
distribution, the sensitivity to 5cp can be improved, and one can discriminate all the values 
of 5cp, including the difference between 5cp = 0 and 7 r . Figure 14 shows the reconstructed 
neutrino energy distributions of the v IJ; sample for several values of 5cp- As expected the 
difference is very small for zv /t events. 

4-4- Analysis method 

The sensitivity of a long baseline experiment using Hyper-K and J-PARC neutrino beam is 
studied using a binned likelihood analysis based on the reconstructed neutrino energy distri¬ 
bution. Both v e appearance and disappearance samples, in both neutrino and antineutrino 
runs, are simultaneously fitted. 

The x 2 used in this study is defined as 

X 2 = —21n£ + P, (14) 

where ln£ is the log likelihood for a Poisson distribution, 

-2In/! = J2 + fi) + Af Ue ln [^(1 + /,:)] } • (15) 

k 

Here, Njf ue (Njf st ) is the number of events in k -th reconstructed energy bin for the true (test) 
oscillation parameters. The index k runs over all reconstructed energy bins for muon and 
electron neutrino samples and for neutrino and anti-neutrino mode running. The parameters 
fi represent systematic uncertainties. For anti-neutrino mode samples, an additional overall 
normalization parameter with 6% prior uncertainty is introduced to account for a possible 
uncertainty in the anti-neutrino interaction, which is less known experimentally in this energy 
region. This additional uncertainty is expected to decrease as we accumulate and analyze 
more anti-neutrino data in T2K, but we conservatively assign the current estimate for this 
study. A normalization weight (1 + /^ orm ) is multiplied to Njf st in the anti-neutrino mode 
samples. 

The penalty term P in Eq. 14 constrains the systematic parameters f t with the normalized 
covariance matrix C, 

P = Yfi{C- l )iJ r (16) 

i,j 

In order to reduce the number of the systematic parameters, several reconstructed energy 
bins that have similar covariance values are merged for f). 

The size of systematic uncertainty is evaluated based on the experience and prospects of 
the T2K experiment, as it provides the most realistic estimate as the baseline. We estimate 
the systematic uncertainties assuming the T2K neutrino beamline and near detectors, taking 
into account improvements expected with future T2K running and analysis improvements. 
For Hyper-K a further reduction of systematic uncertainties will be possible with upgrade of 
beamline and near detectors, improvements in detector calibration and analysis techniques, 
and improved understanding of neutrino interaction with more measurements. In particular, 
as described in Sec. 3.2, studies of near detectors are ongoing with a goal of further reducing 
systematic uncertainties. The sensitivity update is expected in the near future as the near 
detector design studies advance. 
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Fig. 15: Fractional error size for the appearance (left) and the disappearance (right) samples 
in the neutrino mode. Black: total uncertainty, red: the flux and cross-section constrained 
by the near detector, magenta: the near detector non-constrained cross section, blue: the far 
detector error. 


There are three main categories of systematic uncertainties. We assume improvement from 
the current T2K uncertainties for each category as follows. 

i) Flux and cross section uncertainties constrained by the fit to current near detector data: 

These arise from systematics of the near detectors. The understanding of the detector 
will improve in the future, but this category of uncertainties is conservatively assumed 
to stay at the same level as currently estimated. 

ii) Cross section uncertainties that are not constrained by the fit to current near detector data: 

Errors in this category will be reduced as more categories of samples are added to the 
near detector data fit, which constrains the cross section models. We assume the uncer¬ 
tainties arising from different target nucleus between the near and the far detectors 
will become negligible by including the measurement with the water target in the near 
detector. 

iii) Uncertainties on the far detector efficiency and reconstruction modeling: Because most 

of them are estimated by using atmospheric neutrinos as a control sample, errors in 
this category are expected to decrease with more than an order of magnitude larger 
statistics available with Hyper-K than currently used for T2K. Uncertainties arising 
from the energy scale is kept the same because it is not estimated by the atmospheric 
neutrino sample. 

The flux and cross section uncertainties are assumed to be uncorrelated between the neutrino 
and anti-neutrino running, except for the uncertainty of Vejvp, cross section ratio which 
is treated to be anti-correlated considering the theoretical uncertainties studied in [86]. 
Because some of the uncertainties, such as those from the cross section modeling or near 
detector systematics, are expected to be correlated and give more of a constraint, this is 
a conservative assumption. The far detector uncertainty is treated to be fully correlated 
between the neutrino and anti-neutrino running. 

Figures 15 and 16 show the fractional systematic uncertainties for the appearance and 
disappearance reconstructed energy spectra in neutrino and anti-neutrino mode, respectively. 
Black lines represent the prior uncertainties and bin widths of the systematic parameters 
fi, while colored lines show the contribution from each uncertainty source. Figure 17 shows 
the correlation matrix of the systematic uncertainties between the reconstructed neutrino 
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Fig. 16: Fractional error size for the appearance (left) and the disappearance (right) samples 
in the anti-neutrino mode. Black: total uncertainty, red: the flux and cross-section con¬ 
strained by the near detector, magenta: the near detector non-constrained cross section, 
blue: the far detector error. 
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Fig. 17: Correlation matrix between reconstructed energy bins of the four samples due to 
the systematic uncertainties. Bins 1-8, 9-20, 21-28, and 29-40 correspond to the neutrino 
mode single ring e-like, the neutrino mode single ring //-like, the anti-neutrino mode single 
ring e-like, and the anti-neutrino mode single ring //-like samples, respectively. 


energy bins of the four samples. The systematic uncertainties (in %) of the number of 
expected events at the far detector are summarized in Table 9. 

4-5. Expected sensitivity to CP violation 

Figure 18 shows the 90% CL allowed regions on the sin 2 26i^-6cp plane. The results for 
the true values of Scp = (—90°, 0,90°, 180°) are overlaid. The top (bottom) plot shows the 
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Fig. 18: The 90% CL allowed regions in the sin 2 29\3-5cp plane. The results for the true val¬ 
ues of §cp = (—90°, 0, 90°, 180°) are overlaid. Top: normal hierarchy case. Bottom: inverted 
hierarchy case. Red (blue) lines show the result with Hyper-K only (with sin 2 2#i3 constraint 
from reactor experiments). 
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Fig. 19: Expected significance to exclude sin 5qp = 0. Top: normal hierarchy case. Bottom: 
inverted hierarchy case. 
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Table 9: Uncertainties (in %) for the expected number of events at Hyper-K from the 
systematic uncertainties assumed in this study. 


Appearance 

v mode 

Disappearance 


_ Appearance 

v mode 

Disappearance 


Flux & ND-constrained 
cross section 
TO 
2.8 
±6 
4.2 


ND-independent 
cross section 
L2 
1.5 
±0 
1.4 


Far detector 


0.7 

1.0 

L7 

1.1 


Total 


3.3 

3.3 

±2 

4.5 



Integrated beam power [MW 10 7 sec] 


Fig. 20: Fraction of 5cp for which sinJc'p = 0 can be excluded with more than 3 a (red) and 
5 a (blue) significance as a function of the integrated beam power. For the normal hierarchy 
case. The ratio of neutrino and anti-neutrino mode is fixed to 1:3. 


case for the normal (inverted) mass hierarchy. Also shown are the allowed regions when we 
include a constraint from the reactor experiments, sin 2 29^ = 0.100 ±0.005. With reactor 
constraints, although the contour becomes narrower in the direction of sin 2 2 #i 3 , the sensi¬ 
tivity to 5cp does not significantly change because 5cp is constrained by the comparison 
of neutrino and anti-neutrino oscillation probabilities by Hyper-K and not limited by the 
uncertainty of $ 13 . 

Figure 19 shows the expected significance to exclude sindc'p = 0 (the CP conserved case). 
The significance is calculated as \J Ay 2 , where Ay 2 is the difference of y 2 for the trial value 
of 5cp and for 5cp = 0° or 180° (the smaller value of difference is taken). We have also 
studied the case with a reactor constraint, but the result changes only slightly. Figure 20 
shows the fraction of 5cp for which sinJc'p = 0 is excluded with more than 3 <7 and 5cr 
of significance as a function of the integrated beam power. The ratio of integrated beam 
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Fig. 21: Expected 68 % CL uncertainty of 5cp as a function of integrated beam power. 


Table 10: Expected la uncertainty of Aral and sin 2 023 for true sin 2 #23 = 0.45,0.50,0.55. 
Reactor constraint on sin 2 2013 = 0.1 ± 0.005 is imposed. 


True sin 2 023 
Parameter 

0.45 
Am§ 2 (eV 2 ) 

sin 2 023 

0.50 

A mj 2 (eV 2 ) 

sin 2 023 

0.55 
Am§ 2 (eV 2 ) 

sin 2 023 

NH 

1.4 x 10~ 5 

0.006 

1.4 x 10~ 5 

0.015 

1.5 x 10~ 5 

0.009 

IH 

1.5 x 10~ 5 

0.006 

1.4 x 10~ 5 

0.015 

1.5 x 10" 5 

0.009 


power for the neutrino and anti-neutrino mode is fixed to 1:3. The normal mass hierarchy 
is assumed. The results for the inverted hierarchy is almost the same. CP violation in the 
lepton sector can be observed with more than 3(5) a significance for 76(58)% of the possible 
values of 5c p- 

Figure 21 shows the 68% CL uncertainty of 5cp as a function of the integrated beam 
power. With 7.5 MWxl0 7 sec of exposure (1.56xl0 22 protons on target), the value of 5cp 
can be determined to better than 19° for all values of 5cp- 
As the nominal value we use sin 2 023 = 0.5, but the sensitivity to CP violation depends 
on the value of 023- Figure 22 shows the fraction of 5cp for which sin^p = 0 is excluded 
with more than 3 a and 5 a of significance as a function of the true value of sin 2 $23 with the 
current best knowledge of the possible sin 2 $23 range by T2K collaboration [38]. 

4-6. Sensitivity to Am 2 2 and sin 2 $23 

The result shown above is obtained with sin 2 623 and Am 2 2 as free parameters as well as 
sin 2 2013 and 5cp, with a nominal parameters shown in Table 6. The use of the sample 
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Fig. 22: Fraction of 5cp for which sin <5^ = 0 can be excluded with more than 3 cr (red) and 
5 a (blue) significance as a function of the true value of sin 2 023, for the normal hierarchy case. 
Vertical dashed lines indicate 90% confidence intervals of sin 2 023 from the T2K measurement 
in 2014 [38], 


in addition to v e enables us to also precisely measure sin 2 023 and Am| 2 . Figure 23 shows 
the 90% CL allowed regions for the true value of sin 2 023 = 0.5 together with the 90% CL 
contour by T2K v jL disappearance measurement [38]. Hyper-K will be able to provide a 
precise measurement of sin 2 023 and Am| 2 . Figure 24 shows the 90% CL allowed regions on 
the sin 2 023-Am§ 2 plane, for the true values of sin 2 023 = 0.45 and Am| 2 = 2.4 x 10 -3 eV 2 . 
For the determination of 023, a disappearance measurement provides precise measure¬ 
ment of sin 2 2023- However, when 023 % |, there are two possible solutions (023 and | — 023 ) 
which give the same sin 2 2023 • This is known as the octant degeneracy. As seen from 
Eq. 7, u e appearance measurement can determine sin 2 023 sin 2 2013. In addition, the reac¬ 
tor experiments provide almost pure measurement of sin 2 2013. Thus, the combination of 
those complimentary measurements is known to be able to resolve this degeneracy if 023 is 
sufficiently away from ^ [87-89]. As shown in Fig. 24, with a constraint on sin 2 2013 from the 
reactor experiments, Hyper-K measurements can resolve the octant degeneracy and precisely 
determine sin 2 023 - 

The expected precision of Am| 2 and sin 2 023 for true sin 2 023 = 0.45,0.50,0.55 with reactor 
constraint on sin 2 2013 is summarized in Table 10. 

4-7. Combination with atmospheric neutrino data 

Atmospheric neutrinos can provide an independent and complementary information to the 
accelerator beam program on the study of neutrino oscillation. For example, through the 
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Fig. 23: The 90% CL allowed regions in the sin 2 023-Am^ plane. The true values are 
sin 2 023 = 0.5 and Am| 2 = 2.4 x 10 -3 eV 2 (red point). Effect of systematic uncertainties 
is included. The red (blue) line corresponds to the result with Hyper-K alone (with a reactor 
constraint on sin 2 20i3). The dotted line is the 90% CL contour from T2K experiment [38] 
with the best fit values indicated by a black point. 


matter effect inside the Earth, a large statistics sample of atmospheric neutrinos by Hyper-K 
will have a good sensitivity to the mass hierarchy and 023 octant. 

Assuming a 10 year exposure, Hyper-K’s sensitivity to the mass hierarchy and the octant 
of 023 by atmospheric neutrino data are shown in Fig. 25. Depending on the true value of 023 
the sensitivity changes considerably, but for all currently allowed values of this parameter 
the mass hierarchy sensitivity exceeds 3 a independent of the assumed hierarchy. If 023 is 
non-maximal, the atmospheric neutrino data can be used to discriminate the octant at 3 cr 
if sin 2 023 < 0.46 or sin 2 023 > 0.56. 

In the previous sections, the mass hierarchy is assumed to be known prior to the Hyper- 
K measurements. This is a reasonable assumption considering the increased opportunities, 
thanks to a large value of 0 i 3 , of ongoing and proposed projects for mass hierarchy deter¬ 
mination. However, even if the mass hierarchy is unknown before the start of experiment, 
Hyper-K itself will be able to determine it with the atmospheric neutrino measurements. 

Because Hyper-K will observe both accelerator and atmospheric neutrinos with the same 
detector, the physics capability of the project can be enhanced by combining two comple¬ 
mentary measurements. As a demonstration of such a capability, a study has been done by 
simply adding Ay 2 from two measurements, although in a real experiment a more sophisti¬ 
cated analysis is expected. Assuming the true mass hierarchy of normal hierarchy and the 
true value of 5c p = 0, the values of expected Ay 2 as a function of 5c p for each of the 
accelerator and atmospheric neutrino measurements, without assumption of the prior mass 
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Fig. 24: 90% CL allowed regions in the sin 2 023 _ A77?J 2 plane. The true values are sin 2 023 = 
0.45 and Am 2 2 = 2.4 x 10 -3 eV 2 . Effect of systematic uncertainties is included. Top: Hyper- 
K only. Bottom: With a reactor constraint. 


35/40 





8 cp Uncertainty 



sin 2 e 23 



sin % 23 


Fig. 25: Atmospheric neutrino sensitivities for a ten year exposure of Hyper-K assuming the 
mass hierarchy is normal. Top: the Ay 2 discrimination of the wrong hierarchy hypothesis as a 
function of the assumed true value of sin 2 $ 23 - Bottom: the discrimination between the wrong 
octant for each value of sin 2 f? 23 - The uncertainty from 5qp is represented by the thickness of 
the band. Vertical dashed lines indicate 90% confidence intervals of sin 2 $23 from the T2K 
measurement in 2014 [38]. 
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Fig. 26: Combination of the accelerator and atmospheric data. Top: Expected Ay 2 values 
for accelerator and atmospheric neutrino measurements assuming that the mass hierarchy 
is unknown. The true mass hierarchy is normal hierarchy and the true value of 5qp = 0. 
Bottom: By combining the two measurements, the sensitivity can be enhanced. In this 
example study, the Ay 2 is simply added. 
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hierarchy knowledge, are shown in the top plot of Fig. 26. For the accelerator neutrino mea¬ 
surement, there is a second minimum near 5qp = 150° because of a degeneracy with mass 
hierarchy assumptions. On the other hand, the atmospheric neutrino measurement can dis¬ 
criminate the mass hierarchy, but the sensitivity to the CP violating phase 5 cp is worse than 
the accelerator measurement. By adding the information from both measurements, as shown 
in the bottom plot of Fig. 26, the fake solution can be eliminated and a precise measurement 
of 5cp will be possible. 

5. Conclusion 

The sensitivity to leptonic CP asymmetry of a long baseline experiment using a neutrino 
beam directed from J-PARC to the Hyper-Kamiokande detector has been studied based on a 
full simulation of beamline and detector. With an integrated beam power of 7.5 MWx 10 7 sec, 
the value of 5qp can be determined to better than 19° for all values of Scp and CP violation 
in the lepton sector can be observed with more than 3 o (5 a) significance for 76% (58%) of 
the possible values of 5c p- 

Using both v e appearance and disappearance data, a precise measurement of sin 2 623 
will be possible. The expected la uncertainty is 0.015(0.006) for sin 2 623 = 0.5(0.45). 
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